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Abstract 



We present the extension of the Lagrangian loop gauge invariant represen- 
tation in such a way to include matter fields. The partition function of lattice 
<"1 ! compact [/(l)-Higgs model is expressed as a sum over closed as much as open 

surfaces. These surfaces correspond to world sheets of loop-like pure electric 
flux excitations and open electric flux tubes carrying matter fields at their 
ends. This representation is connected by a duality transformation with the 
topological representation of the partition function (in terms of world sheets 
of Nielsen-Olesen strings both closed and open connecting pairs of magnetic 
monopoles). We have simulated numerically the loop action equivalent to the 
Villain form of the action and mapped out the ^-7 phase diagram of this 
model. By virtue of the gauge invariance of this description the equilibrium 
configurations seems to be reached faster than with the ordinary gauge- variant 
descriptions. 
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1 Introduction 



The use of loops for a gauge invariant description of Yang-Mills theories may be 
traced back to the Mandelstam |]J quantization without potentials. In 1974 Yang 
noticed the important role of the holonomies for a complete description of gauge 
theories. 

In the early eighties a loop based Hamiltonian approach to quantum electro- 
magnetism was proposed and generalized [Q in 1986 to include the Yang-Mills the- 
ory. This Hamiltonian formulation was given in terms of the traces of the holonomies 
(the Wilson loops) and their temporal loop derivatives as the fundamental objects. 
They replace the information furnished by the vector potential and the electric field 
operator, respectively. These gauge invariant operators verify a closed algebra (non 
canonical) and may be realized on a linear space of loop dependent functions. 

Afterwards the loop representation was extended in such a way to include matter 
fields: the so-called P-representation This extension includes open paths in 

addition to the closed ones or loops []. 

The loop approach has many appealing features from the methodological point 
of view: First, it allows to do away with the first class constraints of the gauge 
theories (the Gauss law). Second, the formalism only involves gauge invariant ob- 
jects. This makes it specially well suited to study 'white' objects as mesons and 
barions in QCD because the wave function will only depend on the paths associ- 
ated with the physical excitations. Third, all the gauge invariant operators have a 
simple geometrical meaning when realized in the loop space. Additionally, there is 
a conceptual motivation: The introduction by Ashtekar |7j of a new set of variables 
that cast general relativity in the same language as gauge theories allowed to ap- 
ply loop techniques as a natural non-perturbative description of Einstein's theory. 
Furthermore, the loop representation appeared as the most appealing application 
of the loop techniques to this problem ||, @. Thus, it was realized that this loop 
formalism goes beyond a simple gauge invariant description and in fact it provides 
a natural geometrical framework to treat gauge theories and quantum gravity. 

The non-canonical nature of the loop algebra have made elusive a Lagrangian 



loop formalism counterpart of the original Hamiltonian one |H|, pT| . In a previous 



paper |L2| we showed a natural procedure in order to cast loops into the Lagrangian 



formalism. In reference [T2|, it was considered the lattice 4D pure QED model and 



1 Here we will use the term 'loop' for the configurations in presence of matter fields albeit in a 
relaxed sense which covers both closed as much as open paths. 
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it was written the action for the loops as a sum over their closed world sheets. 



In the present paper we shall continue with the program of setting up the La- 
grangian correlative of the Hamiltonian loop formalism. Concretely, we extend this 
Lagrangian loop approach in such a way to include matter fields. 

We consider the lattice compact version of the scalar electrodynamics SQED C 
which describes the interaction of a compact gauge field 0^ with the scalar field 
cf) = \4>\e ltp . The self-interaction of the scalar field is given by the potential A(|0| 2 — 
0o) 2 . For simplicity we shall consider the limit A —>■ oo which freezes the radial 
degree of freedom of the Higgs field (it is known that the numerical results obtained 
already at A = 1 are indistinguishable from the frozen case). Thus the dynamical 
variable is compact, i.e. tp G (— 7r, 7t]. This model is known to possess three phases, 
namely confining, Higgs and Coulomb [I3|. The Higgs phase splits into a region 



where magnetic flux can penetrate in form of vortices (Nielsen-Olesen strings) and 
a region where the magnetic flux is completely expelled ||14|| , the relativistic version 
of Meissner effect in superconductivity. Relying on this, we call this two subregions: 
Higgs I and II in analogy with superconducting materials. 

2 The Loop Representation: Hamiltonian and La- 
grangian Formulations 

The Hamiltonian loop representation of the scalar compact QED in a lattice - with 
sites denoted by x, links denoted by I and plaquettes by p - is given in terms of the 
fundamental gauge operators of the P-representation, the &(P%) defined by 2 

$(PJ0 = &(x)U(I»)4>(y) = ft{x) J] U(l)4>(y) (1) 

(U(l) are the lattice gauge group operators, (j)(x) are the matter field operators) 
and their conjugate momenta: the electric field operator E(l) and the IF (x) and 
tl(x). They obey the commutation relations 



E(i)^(py) =N l (PMpy 



2 Pjf comprises open as much as closed paths or loops C at fixed t, representing "electro- 
mesons" and pure gauge excitations respectively. In this last case $(C) = W{C) i.e. it reduces to 



the Wilson loop operator. 
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n(x),$(p- 



-iS x A(Py) (2) 



where Ni(P) is the number of times the link / appears in the path P. 
The $(-P^) operators are the creation operators of the loops i.e. 

*(J»)|0 >= ft(x)U(P^(y)\0 >= \py > (3) 
where | > is the zero loop state (strong coupling vacuum of the system). 

In order to cast the preceding loop description in the Lagrangian formalism let 
us consider the partition function for the Villain form of the lattice actionP] which is 
given by 

Z = J (dO) E J (dip) E exp(-^ 1 1 W -2vm || 2 -| 1 1 - 2irk - 1 1 2 ) , (4) 



where we use the notations of the calculus of differential forms on the lattice of [15|. 
In the above expression: (3 = \, 6 is a real compact 1-form defined in each link of 
the lattice and ip is a real compact 0-form defined on the sites of the lattice, V is 
the co-boundary operator, n are integer 2-forms defined at the lattice plaquettes, 
and k integer 1-forms, and || . || 2 =< ., • >• 

If we use the Poisson summation formula J2 n f( n ) — J2 S J^°oo d(/)f(4>)e 2m ' l)S for 
each of the integer variables, the partition function transforms into 

(d9) / (dtp) / (#) / (dx) exp(-^ || W - 2^ || 2 

J J—oo J—oo Z 

_| || V ^ - 2n X - 9 ||2) e ^<^> e ^<t, x >_ (5) 
Integrating in the -0 and \ variables 

1 ,,o 1 



z « ee/ m J m^(~ ii s H 2 -- H t ii 2 ) 



Using the partial integration rule < xp, V0 >=< dip, 4> > (d = *V* is the boundary 
operator which maps k-forms into (k-l)-forms and where * is the duality operation 

3 The choice of the Villain form instead of the ordinary Wilson form is only done for simplicity, 
with the Wilson action it is also possible to repeat everything of what follows. 
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which maps k-forms into (4-k)-forms) and integrating over the compact (p and 8 we 
get the constrains 5(dt = 0) and 5(ds = t) and thus, we finally arrive to 



ZcxJ2^M-^\\s\\ 2 -^\\ds\\ 2 ) (7) 



2(3 " " 2k 

or 



i va + M 2 

Z oc > expl < s, — — s > 



2/5 ' M 2 

where M 2 = -| is the mass acquired by the gauge field due to the Higgs mechanism. 
If we consider the intersection of one of the surfaces defined by the integer 2-forms s 
(open and closed surfaces) with ai = constant plane we get spatial paths or 'loops'. 
It is easy to show that the creation operator of those paths is just the the $(PJf) 
operator. Repeating the steps from Eq.(f|) to Eq.flj]) we get for < $(PJf) > 

< HP") E ex P(-^ II * II 2 II 9s + py || 2 ). (9) 

Thus, we have arrived to an expression of the partition function in terms of the 
world sheets of electric string-like configurations: the loop (Lagrangian) represen- 
tation. In this representation the matter fields are naturally introduced by means 
of open surfaces which are the world sheets of open paths representing the 'meson- 
like' configurations. In other words, as it is shown in Fig.l, cutting the above open 
surfaces with planes t = constant we get the corresponding quantum Hamiltonian 
description in terms of open paths Pjf (the P-representation). 



3 Duality and the topological representation 

Another equivalent description of the Villain form is the topological representation 
in terms of the topological objects. As our model has two compact variables we 
have two topological excitations: monopoles and Nielsen-Olesen strings [I6|]. The 
BKT expression for the partition function of compact scalar QED is obtained via 
the Banks-Kogut-Myerson transformation [O (see Appendix) and is given by 



M 2 

Z oc ]T exp(-27r 2 /3 < *n(m), M2 * n(m) >) (10) 

n(m) 
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where m = d * n are closed integer 1-forms attached to links which represent 
monopole loops and *n(m) = *n — d * q are integer 2-forms attached to plaque- 
ttes corresponding to the world sheets of both Dirac and Nielsen- Olesen strings 
(with monopole loops as borders). Thus, comparing (^) and ( ]TD| ) we can observe a 
complete parallelism: in both representations we have a sum over surfaces, and in- 
tersecting with a plane t = constant we get closed as much as open strings with point 
charges at their ends. In the first case this string-like excitations are 'electric' whilst 
in the second they are 'magnetic'. Furthermore, there is a duality transformation 
connecting the confining and Higgs II sectors of the phase diagram. We want to 
remark that there is a slight difference between both equivalent descriptions. In the 
BKT representation monopoles occur at the ends of both the Nielsen-Olesen strings 
(physical excitations) and the Dirac strings (non physical gauge-variant objects) so 
we have the corresponding two types of world sheets mixed in the 2-form *n(m) of 
equation (|T0|). On the other hand the gauge invariant loop description is simpler 
and completely transparent from the geometrical point of view. 

4 Numerical Analysis 

We have performed a standard Metropolis Monte Carlo simulation with the loop 
action of (^). We have worked on a hypercubical lattice. Basic variables are integers 
n^vir) attached to plaquettes p = (r; fi, u), where r denotes the site and (J,, v — 1 — 4 
two directions. These variables describe the world-sheets of the loop-excitations. 

The loop action of (0) is equivalent to the Villain form. There were not numerical 
results of the phase diagram (as far as we know) for the Villain action of SQED C . 
Thus, in order to map out the phase diagram, we made preliminary Monte Carlo 
runs along lines of fixed (3 or 7 on a 6 4 lattice. Typically, 5,000 iterations were 
done at each coupling and measurements of the specific heat were taken. Then the 
coupling was changed by ±0.005, and another 5,000 iterations were made, etc. The 
resulting crude phase diagram (three lines: confining- Coulomb, confining-Higgs and 
Higgs-Coulomb) proved very helpful in guiding our larger scale simulations. 

The next step was to localize accuratelly the transition lines in a 8 4 lattice. With 
this aim, we have performed runs of 40,000 iterations after 15,000 termalization 
sweeps per-point on the plane (3 — 7. The phase diagram is shown in Fig. 2. 

It is worth while to mention that the loop description avoids the problem of sum- 
ming over gauge redundant configurations. This is reflected in the fact that when 
simulating the loop action the convergence to the equilibrium configurations seems 
to be quite faster than when using the ordinary action in terms of the fields. Fur- 
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thermore, we observed the absence of the strong metastability previously observed 
in Monte Carlo analysis of the SQED C lattice theory with the Wilson action, which 
makes very difficult the numerical analysis of the phase transition critical exponents. 
This fact is one of the main advantages for using the action defined in eq. (|?p. 



5 Conclusions 

First, we have shown how to introduce matter fields in the Lagrangian loop rep- 
resentation. It turns out that the partition function of gauge-Higgs theory can be 
represented as the sum over world sheets of loops (open and closed). 

Second, a correspondence between the loop and the BKT descriptions is patent 
and suggests a duality transformation connecting the confining and the Higgs II 
regions of the phase diagram. Perhaps this connection found in the lattice also 
holds in the continuum, supporting the claims on the existence of a new phase in 
QED @. 

Third, simulating numerically a loop-action equivalent to the Villain form of 
the SQED model (with frozen radial degree of freedom) we mapped out its phase 
diagram. It turns out that it is very similar to the corresponding phase diagram 
for the Wilson form [I9|, [20j. Additionally, these results confirm those previously 



obtained for the same model but using the Hamiltonian loop approach |2T| , because 
the mapped phase diagram is qualitatively equal in both cases. 

The next step will be the implementation of the Pauli exclusion principle in the 
Lagrangian loop approach in order to include fermionic fields. This task has been 
accomplished in the context of the Hamiltonian loop formalism in reference || where 
a transparent geometrical description of 'full' QED was given. 



Appendix 

To obtain the monopole representation ([FT]) we start with the Villain form (^) and 
fix the gauge (p = 0. Then, we parameterize the n = Vg + n(v), where q run over 
arbitrary 1-forms and v over all co-closed 3-forms (Vt> = 0). n(v) is a solution of 
Vn = v. If we perform a translation k — > k — q we get the expression 

/+oo R l~ 
(dA) E ex P("- II VA-2ym || 2 || A + 2nk || 2 ) (11) 

-oo i 2 2 

V 

(Vv = 0) 
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where A = 9 — 2irq is a non-compact variable. By shifting A by 2itk we find 
dependence on the combination n + Vk, which turns to be a solution of Vn = v so 
we can eliminate the k variable 

/+oo R K 

{dA) £ ex P(-7T II VA-2vm || 2 -- || A || 2 ) (12) 
-oo * * 2 2 

f 

(V« = 0) 

and performing the gaussian integration we obtain 

M 2 

Zoc eM-^ 2 f3<n(v), vd + M2 n(v) >) (13) 

f 

(Vu = 0) 

Performing a duality transformation we get ( |TUD where m = *v are now integer 
closed 1-forms (dm = 0) and *n(m) = *n(v) + d * q so d * n = m. It is possible to 
express (0) in terms of the k instead of the n variables which reflects the presence 
of Dirac and Nielsen-Olesen sheets. 
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Figure captions 

• Figure 1: Scalar QED surfaces as string propagation. 

• Figure 2: Phase diagram of scalar QED. 



PS files of the figures will be sent by e-mail if 
requested to jomaroca@mat.upc.es 
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